Spin continuity equation and definition of spin current 
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For an electron, a spin-1/2 particle, the spin charge s, a real pseudovector with constant length, 
could determine the spin polarization properties in quantum mechanics. Since spin density ps could 
be expressed as the product of probability density p and the spin charge s, the spin continuity 
equation could be derived from fundamental principles of quantum mechanics and the definition 
of spin current arise naturally. Only in some specific conditions the conventional definition of spin 
current coincides with ours. The equilibrium spin currents would vanish automatically in two- 
dimensional electron gas with Rashba spin-orbit interaction by using our definition of spin current. 
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Spin current is one of the most important physical 
quantities in spintronics. However, there is still not an 
unambiguous definition of spin current [l], 3 • New defini- 
tions of spin current had been proposed [3, 0, S 01 and 
explanations of the conventional definition of spin cur- 
rent had been made 0, B Q ■ Although the interests and 
publications of research in spin dynamics and transport 
still increased in spite of the debates on the definition 
of spin current, the ambiguity of definition of spin cur- 
rent should and must be eliminated as soon as possible. 
In the general transport theory the current density of 
a local charge is defined from the continuity equation. 
It was surprising that spin continuity equation which is 
the continuity equation of spin transport could not be 
easily obtained by computing the time derivative of spin 
density. The reasons come from that the spin continuity 
equation would contain the extra torque terms because of 
the non-conservation of spin and the algebraic properties 
of the spin operator s ~ is complex, where a are the 
Pauli matrices. In quantum mechanics the state of an 
electron, a spin-1/2 particle, is described by the spinor 
*(r, t) and p(r, t) = *1^(r, t)^'(r, t) is the probability den- 
sity. The spin density ps — ^'^(r, i)s5'(r, t) had been well 
defined [lOj. Comparing to the electronic charge density 
Pe = '^^ {r,t)e'^{r,t), one might regard the spin opera- 
tor s as the spin charge which is the local charge in spin 
continuity equation. However, using the general form 
[charge density] = [charge] x [density] , one might intuition- 
ally define the spin charge 



s(r,t) 



,(r,<) ^'t(r,i)s5'(r,t) 



p{r,t) *t(r,t)^'(r,i) 
X^(r,t)sx(r,t) = -n{r,t), 



(1) 



where the n(r,t) — x' {^yt)^x{^i't) is the polariza- 
tion vector that corresponds to the normalized spinor 
X(r,t) =p-i(r,i)4'(r,i) 

In this letter, we show that the spin charge s could 
determine the spin polarization properties of an elec- 
tron. The spin continuity equation would be easily de- 



rived when the spin density is expressed as ps — ps since 
the spin charge s is an real pseudovector with constant 
length. Then the definition of spin current density would 
arise naturally. The equilibrium spin currents in two- 
dimensional (2D) electron gas with Rashba spin-orbit in- 
teraction (SOI) would vanish automatically by using our 
definition of spin current. 

It is well known that the polarization vector n could 
totally determine the density matrix g = ^(1 -I- n ■ <t) 
of an electron in a pure normalized spinor [l3| . The 
polarization vector n was introduced in nuclear magnetic 
resonance (NMR) by Bloch to simplify the treatment of 
spin precession [ij, |l5| . The polarization vector n of an 
electron is a unit real pseudovector in three-dimensional 
space [l6| and could be measured by experiments 17 1. In 



other words, the spin charge s = |n, a real pseudovector 
with constant length, could totally determine the spin 
polarization properties of an electron ^18'|. Moreover, the 
pseudovector s and the spin density ps = ps could be so- 
phistically obtained by using the mathematical language 
Geometric Algebra where the vector s is called spin vec- 
tor [l^, [2^ 21, 22, [2^. We rename s as spin charge 
to emphasize that it functions as the charge in the spin 
continuity equation. 

Since the spin density could be expressed as Ps — ps 
and the spin charge s is a real pseudovector with constant 
length, the spin continuity equation could be easily de- 
rived from the time derivative of the spin density. Below 
we derive the spin continuity equation by an alternative 
method. In quantum mechanics, the fundamental prob- 
ability continuity equation is 

where j is the probability current density. We multiply 
the spin charge s to both sides of Eq. ([2]) 



s(^ 



V-j)-0. 



(3) 



After moving the s into the time derivative and diver- 
gence operators, one could obtain the spin continuity 



2 



equation 



dt 



V • (js) 



P 



ds 

dt 



j ■ Vs, 



(4) 



where the js is a dyadic tensor and the rule of the diver- 
gence of a dyadic tensor is V • (AB) = BV • A + A • VB. 
Our derivation of the spin continuity equation is based 
on the fundamental principles of quantum mechanics. It 
is thus reasonable to define the spin current density 



(5) 



The spin current density is a measurable 2-rank real 
pseudotensor. Our definition of spin current density is 
similar to the definition of spin current density in a mod- 
ified phenomenological Bloch equation 24 1 . 

The two r.h.s. terms in Eq. (|3|) describe the torques. 
The first torque term is just the product of probability 
density p and classical torque ds/dt. The second torque 
term whose physical meaning might not be obvious is the 
dot product of the gradient of the spin charge Vs and the 
probability current density j. Considering an electron in 
a magnetic field H, the Hamiltonian is — ■ H where 
7 is the gyromagnetic ratio. The spin precession satisfies 
the classical equation IJ, |l5 1 



9^ = -^"^" 



(6) 



If the magnetic field H(r) is time-independent, the spin 
charge vector s(r) of the normalized eigenstate is also 
time-independent. The spin charge s(r) is then parallel 
or anti-parallel to the magnetic field H(r) and the torque 
ds/dt vanishes. At that time only the second torque term 
could describe the torque acts on the electron. 

Conventionally, the probability current density j could 
be expressed as j = Re^I'tv^', where v ~ dH/dp and 
H is the Hamiltonian. The z-dircction fc-component spin 
current density is 



(7) 



When the state 5" is the eigenstate of either Vi or Sk , it 
is easily to prove that 



1 



(8) 



which has the same form as the conventional definition 
of spin current density [if. 

It was noticed by Rashba that there would be equilib- 
rium spin currents in the 2D electron gas with Rashba 
SOI by using the conventional definition of spin current 
"if and debates have arisen from then on [1, 0, 0, @, 0j 
@|. Since the state 4" is usually not the eigenstate of 
either Vi or Sk, the i-direction fc-component spin current 
density J'^'^ could not be expressed in the conventional 
form Eq. ([S]). We can show that the equilibrium spin 



currents would vanish automatically by using our defini- 
tion of spin current. In the 2D electron gas with Rashba 
SOI the spin polarization is perpendicular to the k in 
the eigenstates ip),(k), where A = ± correspond to the 
upper and lower branches of the spectrum, respectively 
In any eigenstate the spin charge s is conserved, the 
spin current density is divergenceless and the two torque 
terms vanish. For an arbitrary eigenstate ^p+(k) in upper 
branch, there always exists a corresponding eigenstate in 
lower branch '(/'_(k') which has the same velocity v. The 
spin current densities in the two eigenstates ip+ik.) and 
'0-(k') are opposite to each other since their probabil- 
ity current densities j+ and j_ are equal but their spin 
polarizations s+ and s_ are opposite. It would be easily 
proved that the energy of the two eigenstates (k) and 
'0-(k') are the same. According to the principles of stat- 
ical physics, the two eigenstates V'-i-(k) and '0-(k') have 
the same statistical probability in thermodynamic equi- 
librium. The contribution in total spin current density 
provided by the two spin current density in ip+{\<.) and 
'0_(k') would be zero. Thus in thermodynamic equi- 
librium the total spin current in 2D electron gas with 
Rashba SOI would vanish. 

In conclusion, the spin continuity equation of an elec- 
tron has been derived from the time derivative of the 
spin probability. The spin current density has been de- 
fined naturally from the spin continuity equation. The 
relationship between our definition and the conventional 
definition of spin current is discussed. The equilibrium 
spin currents would vanish automatically in 2D electron 
gas with Rashba SOI by using our definition of spin cur- 
rent. 
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